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Abstract 

We consider the double-elliptic generalisation of dynamical systems 
of Calogero-Toda-Ruijsenaars type using finite-dimensional Mukai- 
Sklyanin algebras. The two-body system, which involves an elliptic 
dependence both on coordinates and momenta, is investigated in de- 
tail and the relation with Nambu dynamics is mentioned. We identify 
the 2D complex manifold associated with the double elliptic system as 
an elliptically fibered rational ("1/2K3 ") surface. Some generalisa- 
tions are suggested which provide the ground for a description of the 
N-body systems. Possible applications to SUSY gauge theories with 
adjoint matter in d = 6 with two compact dimensions are discussed. 



1 Introduction 



Hamiltonian systems of both Calogero-Toda-Ruijsenaars type and spin chains 
can naturally be associated with finite or infinite dimensional groups via 
Hamiltonian reduction. The Hamiltonian reduction procedure produces a 
unified description of such systems in terms of free motion on the group-like 
symplectic manifold restricted to a fixed level of the moment map [p]]. The 
line of generalisation from rational to trigonometric to elliptic models has 
a clear algebraic counterpart in terms of zero loop, one-loop and two-loop 
structures M. One further has to distinguish two different possible generalisa- 
tions, where either coordinates or momenta take values on the elliptic curve. 
It is useful to introduce the notion of "dual" system here || |J, where the 
coordinates and action variables effectively get interchanged. The attempt to 
give a mathematically rigorous definition to the duality of integrable systems 
has been recently given in ||. 

The line of generalisation in terms of increasing the number of loops to 
produce say the trigonometric Ruijsenaars models from the rational ones is 
not the only one that has been considered. The same results can be obtained 
for the rational and trigonometric Ruijsenaars models moving in another 
direction ||. Namely, there it was shown that these models admit the reali- 
sation as motion on a finite dimensional quantum group. However a similar 
interpretation for the elliptic models was lacking. Below we shall suggest 
a realisation of the elliptic models as motion associated with certain finite 
dimensional algebras.Q It seems that the proper algebraic object for the gen- 
eral double-elliptic system is the generalised Mukai-Sklyanin algebra. This 
was introduced for two quadrics in C 4 in ]7| and may be easily generalised to 
n quadrics in C n+2 in such a way that the famous Mukai symplectic Poisson 
brackets are obtained by considering the corresponding K3 surfaces as affine 
symplectic leaves of the n quadric structure || ||. 

In this paper we shall use the algebra induced by four quadrics in C 6 which 
yields the phase space for the two-body double-elliptic problem. This model 
was given explicitly in Q in Darboux variables. We shall show that this 
formulation of the double elliptic system is equivalent to the one suggested 
in Q. All other two-body elliptic models such as the elliptic Ruijsenaars 
and Calogero models, as well as the periodic Toda chain, can be derived 
from the double-elliptic system via one or other degeneration. Recall that 
the elliptic Calogero system is an example of a Hitchin system [ITS]. More 



general integrable systems corresponding to Hitchin systems for curves with 
many marked points were discussed in [|14|, [13|, [15| . 



Of course the real problem to be solved is a derivation of the integrable 
many-body system. In the relatively simple cases known so far, one could 
just consider a higher rank group, which can be finite or infinite dimensional. 
This program works perfectly up to the elliptic Ruijsenaars model which is 
the most complicated system treated successfully in this way |TJ|. However 
the next step to the double-elliptic system appears to be very complicated 
technically and so one has to look for another technique. One method, which 



lr The preliminary results of this paper were presented in the review ||. 
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suggests a complementary description, is the separation of variables. This 
approach yields a parametrization of the phase space of the many-body sys- 
tem as the symmetric power of some complicated two dimensional complex 
manifold IJT7J. The most difficult part of analysis in this picture is the proper 
identification of the two-dimensional complex manifold corresponding to the 
two-particle problem. 

The problem of a geometrical description of the double-elliptic systems, 
as well as being intrinsically interesting from the point of integrability, has 
important applications to the description of low-energy effective actions in 
SUSY gauge theories in different dimensions. Via Seiberg-Witten theory it is 
known that the classical dynamics of integrable many-body systems addresses 
quantum issues in SUSY gauge theories (see jl8j [19|] and the references therein 



for a review). The most complicated SUSY systems treated in this manner 
thus far have been the elliptic Ruijsenaars model, corresponding to a D=5 
SUSY theory with adjoint matter and with one compact dimensionpO], |21~H , 



and theories with fundamental matter in D=5, 6 [22 1. The double-elliptic 
system supplies the next step to a D=6 SUSY system with adjoint matter, 
and where two dimensions are compact. The two elliptic moduli to be found 
in the gauge theory can be identified with the complexified coupling constant 
and the ratio of the radii in the fifth and sixth dimensions. Therefore the 
answers we obtain will be of the immediate use for this theory. The pure D=6 
SUSY theory with two compact dimensions corresponds to the degeneration 
of the double-elliptic system to the "elliptic Toda" one. 

Let us also note that a clear representation of the integrable systems cor- 
responding to D=6 SUSY gauge theories is also necessary to formulate the 
dualities from || |]] precisely. The point is that these dualities relate gauge 
theories in different dimensions. For instance, the systems dual to the elliptic 
Calogero or periodic Toda system (D=4) correspond to D=6 gauge models. 
Another example involves the relationship between Ruijsenaars (D=5) mod- 
els and the trigonometric Calogero (D=4) ones. It seems that this duality 
will be very helpful in the nonperturbative (De)-construction of dimensions 
recently discussed in [f23 . 



The paper is organised as follows. We start with general remarks con- 
cerning the generalised Mukai-Sklyanin algebra. Then we identify the four- 
dimensional manifolds responsible for the elliptic models including the double 
elliptic one. The possible generalisation to the n-body problems will be sug- 
gested. Finally we apply the double-elliptic system to the description of 
the SU{2) SUSY gauge models with the adjoint matter in D=6 where two 
dimensions are assumed to be compact. 

2 Two-body system 
2.1 General facts 

The generalised Mukai-Sklyanin algebra provides a Poisson structure on the 
intersection of n polynomials Qi in C n+2 . The case of n = 2 quadrics cor- 
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responds to the Sklyanin algebra while the case n = 3 corresponds to the 
Poisson structure attributed to a K3 surface. The Poisson bracket in affine 
coordinates looks as follows: 



dQ 1 8Q 



(1) 



dx kl ' ' ' dx kn ' 



where affine coordinates in C n+2 . The corresponding four dimen- 

sional manifold is "Poisson noncommutative" , a point we shall amplify in 
due course. The polynomials Qi themselves yield the Casimirs of the alge- 
bras. 

We choose the following system of four quadrics in C 6 which provides the 
phase space for the two-body double-elliptic system 



Here the first pair of equations yields the "affinization" of the projective 
embedding of the elliptic curve into CP 3 , while the second pair provides 
the elliptic curve which locally is fibered over the first elliptic curve. If 
the coupling constant g vanishes the system is just two copies of elliptic 
curves embedded in C 3 x C 3 . Let us emphasise that the coupling constant 
amounts to an additional noncommutativity between the coordinates beyond 
the standard noncommutativity of coordinates and momenta. 

We should also warn that we will use the same letters to denote the 
homogeneous and non-homogeneous coordinates in all non-confusing cases. 

2.2 Poisson structure 

The relevant Poisson brackets for this particular system of quadrics reads 

{x 1 ,x 2 } = {x 1 ,x 3 } = {x 2 ,x 3 } = 0, 

{x 5 , Xi} = -X 2 X 3 X 4 Xq, {X bj X 2 } = -£1X3X4X6, {x 5 , x 3 } = -XiX 2 X 4 X 6 , 

{x 5 , x 4 } = -g 2 Ar 2 2xix 2 x 3 x 6 , {x 5 , x 6 } = 0. 

(Various factors of two have been incorporated into the definitions of the 
quadrics to avoid unnecessary factors in these Poisson brackets.) We should 
remark that the Poisson structure is singular and can't be extended up to 
a holomorphic structure on the whole CP 6 because of the arguments of 0. 
Namely, let Xi, ...,X n be coordinates on C n considered as an affine part of 
the corresponding projective space CP n with the homogeneous coordinates 
(xo : Xi : ■ • • : x n ) , = Then if {Xi, Xj} extends to a holomorphic Pois- 
son structure on CP n the maximal degree of the structure (= the length of 
monomials in Xi) is 3 and 



X fe {X i? X 3 } 3 + X t {X v X k } 3 + X,{X k , X t } 3 = 0, i ^ 3 ^k. (3) 
Thus {Xi,Xj} 3 = XiYj - XjY i7 with degYi = 2. 



1 2 — 

r~f> 2 ,-y»2 

1 3 — 

— g 2 x\ + x\ + x\ = 1 
— g 2 x\ + x\ + k~ 2 x\ = k~ 2 . 



(2) 
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The nontrivial commutation relations between coordinates on the distinct 
tori correspond to the standard phase space Poisson brackets while the non- 
trivial bracket {X5, £4} is the additional noncommutativity of the momentum 
space mentioned earlier. 



2.3 Equations of motion 

Let us note that the triple x±,X2, x 3 can be considered in the elliptic parametriza- 
tion 



^"1 sn(q\k) 

_ cn(q\k) fr\ 

_ dn(q\k) /gN 

•J sn(a\k\ \ ' 



with the other three coordinates of (0) expressed by the natural Jacobi func- 
tions uniformizing the second ("momenta") elliptic curve: 

x 4 = asn(J3p\f) (7) 

x 5 = acn((3p\f) (8) 

XB = /3dn(J3p\f). (9) 

Here 

«M = <l 1 + dwr < 10) 




The Hamiltonian of the double-elliptic system in the form of |4]] 

H{p, q) = a(q\k) cn(pP(q\k, *)|^#) (12) 

coincides^] with x 5 . 

We can check directly (see the Appendix) that the canonical equations of 
motion can be re-written as the following polynomial system 

Xi = {xi,X 5 } = X2X 3 X 4 X 6 

£2 = {^2,^5} = x 1 x 3 x i x e 

x 3 = {x 3 ,x 5 } = xix 2 x 4 a; 6 , , 

x 4 = {x 4 ,x 5 } = g 2 x 1 x 2 x 3 x 6 

%5 = {^5,^5} = 
x s = {x 6 ,x 5 } = 0. 

2 We have chosen a value of the coupling constant g to be pure imaginary and rescaled 
by y/2 with respect to the choice of H . 
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Due to the commutation relations x§ is a constant, Xq = c say, which has 
to be related to the energy x$ = E to provide the consistency of the system 
which comes from the last two equations of (0): 

~k 2 (l-E 2 ) = 1-c 2 . (14) 

This consistency condition reduces the number of quadrics in the system 
because of the coincidence of the last two equations. 

To get an explicit form of the spectral curve which corresponds to the 
solution to the equations of motion let us put y = Xix 2 x 3 x 4 and obtain the 
following equation 

y 2 = x\{x\ - l){x\ - k 2 ){g 2 x\ - (E 2 - 1)). (15) 

Taking s = x\ and assuming the energy level to be "generic" [E 2 ^ 1, g 2 + 
l,g 2 k 2 + 1) we obtain that the curve under consideration is a two-sheeted 
covering of the following hyperelliptic form of an elliptic curve 

y 2 = s{s _ 1){s _ k 2 ){g 2 s _ {E 2_ 1)) 

with the branching points 0, 1, k 2 , (E 2 — l)/g 2 - We can immediately integrate 
the system by the hyperelliptic integral: 

2x\dx\ 2x\dx\ ds 
dt = — - = — — - = — , 

X1X2X3X4 y y 

ds ^ f dxi 



V J ^{x 2 - l){x 2 - k 2 )(g 2 x 2 - (£ 2 - 1)) ' 

Remark 1. Let us suppose that the coupling constant g in the hamilto- 
nian (|T2"D is equal to (the case of "free" motion). Then the expression of 
fli~2f) reduces to 

H(p,q) = cn(p\k) (16) 
and the parametrization of the second curve is simplified as 

£4 = sn{p\k) (17) 
x 5 = cn(p\k) (18) 
Xq = dn{p\k). (19) 

We have the following "decoupled" family of the quadrics in C 3 x C 3 : 

, 2 , yi 2 1 

1 2 — 

, yi 2 ,-vt 2 7^,2 

1 3 — 

x\ + x\ = 1 
k 2 x\ + x\ — 1. 

The geometrical picture of this family is of course an affine part of the 
direct product of the elliptic curves £ q {k) parametrized by x 1; x%, x 3 and £ p {k) 
parametrized by x 4 , x 5 , x 6 . 
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The hamiltonian again is identified with x§ and we obtain a new involutive 
coordinate X4. The canonical equations of the "free" motion may be written 
down immediately: 

X\ = {xi,X$} = X2X3X4X6 

X2 = {^2,^5} = xiXsx^xq 
x 3 = {x 3 ,x B } = x 1 x 2 x i x G 

X4 = {x 4 , £5} = 

£5 = {^5, £5} = 
x 6 = {x 6 ,x 5 } = 0. 

and after restriction on a level 25 = E we can easily recognise in the above 
system the classical Arnold- Euler-Nahm top for SU(2) in the variables x\, X2, x 3 . 

Remark 2. Let us consider the "singular energy level" E = —1. Then 
we have the following nodal rational curve 

y 2 = g 2 s 2 {s-l){s-k 2 ) 

which may be reduced to the following canonical form 

Y 2 = g 2 (l- S)(l-k 2 S) 

by the transformation 

S=l/s, Y = y/s 2 . 

Remark 3. We note that the spectral curve of the double elliptic system 
(|H)|) has a clear physical meaning from the point of view of D=6 gauge 
theory. It is this curve that provides the information concerning the low 
energy effective action as well as the spectrum of BPS states. The energy 
level E corresponds to the coordinate on the Coulumb branch of the moduli 
space in the corresponding gauge theory. The masses of the BPS particles 
can be derived calculating the integrals of the meromorphic differentials over 
the cycles on the spectral curve. The singular values of the energy, at which 
points the curve degenerates, correspond to the singularities on the Coulumb 
branch of the moduli space when the mass of some BPS particle vanishes. 



2.4 Another choice of the double-elliptic ansatz 

The choice of another solution for the ansatz of |4|] should give an another 
reparametrisation for the double-elliptic phase space. It is clear that all such 
choices are related to the initial by an appropriate modular transformation 
of the Jacobi moduli k,k,k e ff,k e ff (see ||). In fact, we are able to check 
the following statement: The second choice of the double-elliptic ansatz in 
0] provides us with a polynomial Hamilton system which is equivalent to the 
first. 

Indeed, the Hamiltonian of the double-elliptic system in the form of the 
second solution of M is 

H(p, q) = a(q\k) cn(p\ka(q\k)) (20) 



7 



where now a(q\k) = w 1 + ^i^y Taking the the following parametrizations 
of the coordinate and momentum tori 

x = 1 x = s "( g l fc ) x = dn ^ k ^ ■ 

1 cn(q\k)' 1 ^ cn(q\k) ' ^ cn(q\k) ' 

y x = sn(p|fca(g, fc)), y 2 = cn(p|fca(g, fc)), y 3 = dn(p\ka(q, k)) 
and then denoting by 



£4 = 2/3, £5 = a/1 + # 2 x 2 y x , x 6 = a/1 + c/ 2 xf y 2 
we obtain the four quadrics in C 6 : 



/ y> 2 /y» 2 "1 

/ >T 2 ,-VT 2 /fS/yiS 

■ jj -y 3 — 2 

— g 2 x\ + %1 + %1 = l 

2^2^,2 _i_ ™ 2 — 1 
rt- 1.0 ^ — J- • 

Again the double-elliptic hamiltonian coincides with 25 and the system 
associated with the quadrics has a form similar to (|13|): 

i?i = {xi,x 5 } = -22X324X6 

X 2 = {£2,0:5} = -£1X324X6 

£3 = {^3, £5} = —k' 2 xix 2 x^x 6 
x 6 = {x 6 ,x 5 } = -g 2 x 1 x 2 x 3 x 4 
x 5 = {x 5 ,x 5 } = 

£4 = {^4, £5} = 0. 

Here k' 2 + k 2 = 1 is the complementary modulus and the variable 24 plays 
the same role as xq in the case (|T3|) and also can be eliminated to obtain 
an almost identical polynomial form. As was argued in this choice of 
solution better the study of degenerations of the model (see below chapter 
4). 



3 Duality 

Let us now describe the "duality" property for the dynamical system under 
consideration. We shall formulate the duality as an involution of the four 
dimensional manifold. First perform the following involution of homogeneous 
coordinates 

xi^x , Xq^-xi, x 2 — > £5 (21) 
which yields the following intersection of four quadrics 

222 

rr> rr> rr> 

5 — 1 

, >> 2 ,->-»2 — 7^2 ,y,2 

— 5» 2 Xq + x\ + x\ — x\ = 

z~ 1 2 2 ™ 2 1 7^2^y,2 1 ™ 2 — ,-y»2 

The dual Hamiltonian of j|] is 

= cn(#) = — (22) 
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In terms if the affine coordinates 1 < % < 6: 

Xo X5 X3 X4. X2 Xq 

Z\ = , Z2 = , Z3 = , Z4 — , Z5 = , Zq = . 

X\ X\ X\ X\ X\ X\ 

our "dual" system of quadrics has "affine" form 

z 2 -z 2 = l 

7 2 _ 2 _ p 

-g 2 z\ + z\ + 4 = 1 
-g 2 k 2 z\ + k 2 z\ + z\ = 1. 

The dual Hamiltonian is encoded in the variable Z2 and corresponding 
Mukai-Sklyanin Poisson algebra gives the "dual" Hamilton system 

Z\ = {Zi,Z 2 } = —Z3Z4Z5ZQ 

z 2 = {z 2 ,z 2 } = 

z 3 = { z 3, z 2 } — —Z1Z4Z5ZQ 

Z4 = {z4,z 2 } = —g 2 ZiZ^z^ZQ 
Z5 — {^5,2:2} = —Z1Z3Z4Z6 
z 6 = {z 6 ,z 2 } = 

which is clearly equivalent ("self-dual") to our initial double-elliptic system 



Remark It is worth noting that the dual Hamiltonian ( p2|) in the initial 
("canonical") coordinates after the Fourier rotation p — > — q, q — > p goes to 
the Hamiltonian ( |T6D of the "free" (g = 0) model corresponding to an Arnold- 
Euler-Nahm SU(2) top. Below we will show that there is a transformation of 
the "non-free" hamiltonian system with general double-elliptic Hamiltonian 
(12) into a system describing a pair of identical SU (2)-Arnold-Euler-Nahm 
tops. 



3.1 Algebro-geometric remarks 

The system of four quadrics (H) can be considered as an affine part, S s , of 
the intersection surface S s C CP 6 : 

222 

X\ X 2 — Xq 

rp 2 ,-y>2 — ^2^2 

1 3 — 

fl 2 ,-yi 2 I ,->t2 I ,y.2 ™ 2 

fi 2 ,-y, 2 I ™ 2 I iU 2 ™2 — h, 2^,2 

y j/ j ~ r «x 4 f J- q — Ay 0* 

This is a singular complex surface (i.e. a two-dimensional algebraic variety) 
which depends on the vector s = (1, k 2 ,g 2 , k~ 2 ). A straight forward compu- 
tation shows that (up to a non-zero factor) the Jacobian for the above system 
is given by the following matrix: 



/ 


x x 


—x 2 














-x \ 




X\ 





-X3 











— k 2 x 




-g 2 xi 








X4 


x 5 







V 














-25 


k~ 2 x 6 


(k- 2 - l)x J 
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There are therefore the two following sets of singular points of S s : 



1) xq = X5 = xq = (J; x 4 = g Xi, x 1 = x 2 = x 



3) 

which contains 8 points and 

2) x = xi = x 2 = x 3 = 0; x\ + x\ = 0; x\ + k~ 2 x\ = 0, 

containing another 4 points. 

We will argue below that our intersection surface belongs to the class 
of so-called Adler-van Moerbeke surfaces in CP 6 which project to a (singu- 
lar) Kummer surface, which in its turn gives a Del Pezzo (elliptically fibered 
rational) surface with 12 singular elliptic fibers ("1/2K3" surface). Our phi- 
losophy here is to "embed" the 1- dimensional integrable system, associated 
with a complex curve (in our case this is the "base" elliptic curve of the co- 
ordinates Xi, X2, X3) into a system associated with a 2-dimensional complex 
variety along with the deformation of Hitchin systems proposed in [|34] . 

Our model has evidently only one degree of freedom so we are interested 
in a section of the surface given by the "energy level" H = z, which is a 
hyperplane section in our case. The canonical divisor K$ = (4x2— 4 — 3)H = 
H is given by a hyperplane section H of the surface S s . 

We will be interested in the hyperplane section C z of the surface by x 5 = 
zx , given by a fixed ("generic") complex number z : C z = S s f]{x 5 = zxq] C 
CP 5 . 

Let us compute the genus of this (singular) curve taking into account 
that it passes through the 8 of the twelve singularities enlisted above. The 
canonical divisor is calculated as = (8 — 6)H = 0^ z {2) and degC 2 = 
{{CD + (C z , Ks) = 2 5 . We have by the adjunction formula 

g(C z ) = 1/2((C 2 2 ) + (C z , K s )) + 1 - 5 = 2 4 + 1 - 8 = 9 

and we have that the genus of the (normalised) curve is 9. 

Let consider the corresponding ("projected" )family of curves in CP 4 
because the last equation of the quadrics is a superficial one due to the 
consistency condition between x 5 and x e (see 2.3). 

The family C z represents a singular family of plane curves 

C z : y 2 = x 2 (x 2 -l)(x 2 -k 2 )(g 2 x 2 -(z 2 -l)), 

with x = x\, y = X1X2X3X4. This family is incomplete and has a double point 
at the origin for each "generic" value of z. A completion of the "generic" 
curves in the family C z may be achieved in two different ways which we will 
now briefly discuss. 

The first is by adding a point at infinity which is a completion of the 
affine curve in C 2 to a singular curve C z C CP 2 of genus g(C z ) = 2. We have 

C z : ^^^(X 2 -^ 2 )^ 2 -^ 2 )^ 2 -^ 2 -!)^ 2 ), 
with (X : Y : Z) homogeneous coordinates on CP 2 such that x = y,y = 
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The curve C z is a double covering of CP 1 , ttq '■ C z — > CP 1 , branching at 
6 points. Its normalization E z has genus 1 because the family C z is a double 
cover of a family of elliptic curves (an "elliptic pencil"): 

E z :Y 2 = (l- X)(l - k 2 X)(g 2 - (z 2 - 1)X), 

which is in turn a double cover of CP 1 , tie '■ E z — > CP 1 with the invariants 

g 2 (z) = 1/12{(A; 4 - k 2 + \){z 2 - l) 2 - k 2 g\k 2 + l){z 2 - 1) + k 4 g 4 } G O z (4); 

g 3 (z) = l/48k 2 (z 2 - l)[{z 2 - 1) + (k 2 + l)g 2 )[(k 2 + l){z 2 - 1) + fcV]- 

-1/216[(A; 2 + 1)(^ 2 - 1) + A;V] 3 - g 2 /\Qk\z 2 - l) 2 G £> 2 (6). 

For the second we can consider the complete curve C z of genus 3 which 
is a normalization of the fiber product 

C z x CP i E z := {(p,q) G C z x P 2 |tt c (p) = 7r E (g)}. 

The curve C z is an "etale double covering" of C z . 

Let us discuss the requirement for the intersection surface S(S S ) to be an 
abelian surface. In what follows we will suppose that it is. The curve C z 
is irreducible of genus 2 curve with one double point admitting an elliptic 
involution. It lies on the surface S. We describe the beautiful geometric 
relations between the triple of the families of curves C z , C z and E z and the 
surface S which is based on pi}] . 

Fix the double covering ir : C z — > E z and let be a branch point. Then 
we have the natural mapping C z — > J{C Z ) := Pic°(C z ) by p G C z i— > Oq z (j> — 
a ). If 7r(a ) = qo then the elliptic curve E z is identified with Pic°(E z ) by the 
mapping q G E z f— > OeXi ~ Qo) so we have the induced mapping 7r* : P 2 — > 
J{C Z ) such that ker7r* = 0. The compatibility of the involutions shows that 
the images of Pic°(E z ) and Pic°(C z ) in Pic°(C z ) are complementary. The 
quotient of Pic°(C z ) £ J(C' Z ) by the subgroup ker(P2 C °(P z ) -> Pic°(C z )) 
is an abelian surface which is called the Prym variety of the involution on 
C z and denoted by Prym(C z /E z ). Then the surface S is identified with 
the quotient J(C z )/ir*(E z ) and is dual to Prym(C z / E z ) by the "duality 
theorem" (th.1.12 in |24|j ). Moreover, both surfaces S and Prym(C z / E z ) are 
(1,2)- polarized and the following diagram takes place: 





Prym{C z /E z ) 




CP 1 
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Here the symbol ~ means that the surfaces Prym(C z / E z ) and J(C Z ) are 
isogenous. The fact that Prym(C z / E z ) ~ J{C Z ) is rather standard (in the 
framework of integrable systems we learnt it from the beautiful book [[37]). 
The symbol •<=>- above means the duality between the intersection surface 
S and the Prym variety Prym(C z / E z ), namely S v -<=>- Prym(C z / E z ) and 
(Prym{C z / E z )) y = S with the duality of (l,2)-polarizations. 

On other hand, it is well-known (see |[25|| ) that the Jacobian of a genus 
2 curve of "type I" (which is here the case for C z which has affine equation 
y 2 = (x 2 — l)(x 2 — k 2 )(g 2 x 2 — [z 2 — 1))) is isogenous to a product of two 
elliptic curves. Hence the dual of the intersection surface S is isogenous to a 
product of two elliptic curves. Naively this isogeny can be interpreted as the 
combination of the elementary duality transformation of |4| and the Fourier- 
Legendre rotation applied to the surface S given by the quadric intersection. 

Remark The "non-generic" values of z (z 2 = 1, g 2 + l,g 2 k 2 + 1) reduce 
the family C z to different examples of plane nodal rational curves. 

3.2 Comparison with Beauville-Mukai systems 

When we fix the consistency condition 

k 2 {l-E 2 ) = 1-c 2 (23) 
a system of the following form emerges: 

/y»2 ,-y>2 1 

1 2 — 

, yi 2 , yi 2 I-.2 

1 3 — 

— g 2 x\ + x\ + x\ = 1. 

We would like to associate the intersection surface S(S S ) with a double 
covering tv : S *—>■ CP 2 a la Beauville-Takasaki ( p6| , p7[). 

Let us suppose (according to the discussion in 0) that the value of the 
Hamiltonian x$ = z is a coordinate function on the open set Uq = CP 1 \ oo 
and take homogeneous coordinates (X : X 1 : X 2 ) on CP 2 . Consider the 
surface S C CP 2 x A 1 (A 1 is an affine part of CP ) given by the equation: 

X X 2 = Xl + f(z)X 2 X 1 + g{z)Xl 

where f(z) and g(z) are polynomials. Again let us put y = x\XiXzX± and 
obtain the equation 

y 2 = x\{x\ - l)(xl - k 2 ){g 2 x\ - (z 2 - 1)). (24) 

from our system of quadrics. With help of the evident linear change 

we can replace the curves of the cubic family 

Y 2 = (1 - X)(l - k 2 X)(g 2 - (z 2 - 1)X). 
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by an equation of Weierstrass form. If we consider the coordinates Y, X as 
non-homogeneous for (X : Xi : X 2 ) (ie Y = X 2 /X , X = Xi/Xq) the 
equation acquires the Weierstrass form: 

Y 2 = X 3 + f(z)X + g(z), (25) 

where f(z) E 0(A),g(z) E 0(6). 

Let us produce the transition to the chart U\ = CP 1 \ with coordinate 
1/z: 

z=l/z,X = X/z 2 ,Y = Y/z 3 . 

We have 

Y 2 = X 3 +~f(~z)X+~g(~z) 1 (26) 

where f(z) = f(z)z~ 4 ,g(z) = g(z)z~ 6 . We can consider this as a surface 
Si C CP 2 x A 1 (where A 1 = Ui) and hence the union So [j Si correctly 
defines a surface / : S E CP 1 which is elliptically fibered. In general it is 
a singular surface (at the points where the J— invariant J(z) = f(z) 3 /A(z), 
with A = Af(z) 3 + 27g(z) 2 E 0(12) vanishes). 

Therefore, we can suppose that the family of curves C z covers the elliptic 
pencil lying in (an open part of) a rational elliptic ("1/2 K3") surface S. 
The rational elliptic fibrations (del Pezzo surfaces) can (roughly speaking) 
be considered as a Kummer surfaces /C modulo Z2- The Kummer surface in 
our case is exactly the 2:1 projection of the intersection (|3.1|): 



c z 


cS s 


c CP 6 


i 


i 


I 


c z 


c K 


c CP 5 


I 


I 




E z 


C K/Za 





Remark This description fits well to the results of [ I I | where the full 
description of Kummer surfaces in CP 5 arising as projections of the Adler- 
van Moerbeke [[12] family of abelian surfaces in CP 6 . Our description is 
corresponded to the first type of the surfaces : a Kummer surface with 12 
nodes. The explicit quadratic equations of all 9 types of singular Kummer 
surfaces with 1:1 correspondence to the Adler-van Moerbeke surfaces can be 
found in []11]J. 



Let us give another geometric picture of the intersection surface. In 
"physical" language we are dealing with a family of "supersymmetric cy- 
cles" (C Z ,L) in S, where C 2 := [C Z ].[C Z ] = 2g(C z ) -2 = 2 (because the 
genus g{C z ) = 2), [C z ] its homology class in H 2 (S,Z) and L is a flat U(l)- 
bundle (:= a choice of a point on J{C Z )). We consider the linear system 
\C Z \ := F{H°{S, 0{C Z )) = CP 2 of all such curves. 

The family C z is singular and we consider a compactification J{C Z ) of the 
Jacobian variety J(C Z ) of C z . This compactified Jacobian is a fibration 

J{C Z ) - J{E Z ) 
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with singular rational fiber. Here J(E Z ) is the Jacobian of the normalisation 
such that the generalised Jacobian is the extension 



-> C* -> J(C Z ) -> J(^) -> 0. 

The union [J J(C Z ) has dimension 2g(C z ) = 4 and is a fibration over |C Z | = 
CP 2 . 

In our case 0{C Z ) generates Pic(S) so the family of the compactified 
Jacobians carries a hyperkahler structure, moreover it is birationally sym- 
plectomorphic to Hilb™(S) - the Hilbert scheme of 2 points on the surface 
S (see [PE| , P^ , [T7[]). This is a framework of the so-called Beauville-Mukai 
integrable systems. Namely we have a holomorphic Lagrangian fibration 
U J(C Z ) — ► |C Z | and an open embedding of Symm 2 (S) — > (J J{C Z ) (choosing 
a base point on C 2 we can use the Abel-Jacobi map AbJ : Symm 2 (C z ) — > 
J(C Z ), AbJ(pi,p2) = Oc z {pi + T>2 — 2poo)), such that the following commu- 
tative diagram takes place 

Symm\S) >-{JJ(C z ). 




Hamiltonian in our case implies that we have a degree 6 polynomial defin- 
ing a hyperelliptic curve of genus 2: 

Y 2 = (cz 2 + u x z + u 2 f + f{z){cz 2 + mz + u 2 ) + g{z). 

Here we should consider the constant c as a Casimir following Takasaki's 



prescription in |27| . Various expressions of the canonical form are considered 
in Appendix B. 

3.3 Comparison of two descriptions of the double-elliptic 

system 

Let us now compare the above description of the double-elliptic system with 
that suggested in ||. The two dimensional manifold which provides the phase 
space in that paper was identified with the double covering of an elliptically 
fibered two-dimensional manifold (K3 surface). The suggested Hamiltonian 
was the coordinate on the base of the fibration while the Hamiltonian of the 
dual system was identified with the coordinate of the fiber. It was supposed 
in H that generically the K3 surface contained an incomplete hyperelliptic 
genus five curve, and that is why Takasaki [p7|] in order to obtain an embed- 



ding of the one degree of freedom double-elliptic system in to the scheme of 
Beauville had to set four of his five parameters itj, i — 1, . . . , 5 to be equal 
zero. But, as was shown above, the manifold of [§] has the structure of a 
double covering of an elliptically fibered rational surface. Moreover in these 
terms the Hamiltonian of the system is the coordinate on the fiber, while the 
Hamiltonian for the dual system is the coordinate on the base. Therefore 
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we could consider the same argument as for the double-elliptic system from 
H in the case of the rational elliptic fibrations, (incomplete) hyperelliptic 
curves of genus 2 (a typical fiber of the dual system) and the double cover 
looks like a hyperelliptic Jacobian of the genus 2 curve. This gives us a direct 
embedding of our picture in the Beauville-Mukai systems a la Takasaki p7[ 



taking the Hilbert scheme (or rather Symm 2 (S)) with the symplectic 
form teWfdzW + dxmjdzm and the Abel _ Jacobi map 

/ r pi zdz r P2 zdz r pi dz r P2 dz\ 

AbJ( Pl , P2 )=( Z -^+ ™; ^+ -)=(0 1 ,0 2 ) 

\jpoo y j Poo y J poo y Jpoo v / 

and 

^ T «{dxWAdzW dx^Adz^\ , 1± , 
AbJ 7-t 1 TTyv = dui A dtpi + du2 A d<p2- 

\ y {1) y [2) ) 

Here we suppose that pi = (z^ l \y^) G S, i = 1, 2 We can illustrate the situa- 
tion by the following diagram of double-elliptic system with the interchange 
of the direct and dual systems. 



S ~ J(C X ) J(C Z ) ~ S v 

x e CP 1 x g CP 1 z e CP 1 z e CP 1 

Moreover, now we are able to explain the "naive" definition of duality 0] 
arising from the anticanonical condition: 

dP AdQ = -dp A dq. 

Consider two "dual" elliptic fibrations associated with the double-elliptic 
system and with its dual in terms of the families of projective Weierstrass 
cubics in CP 2 : 

S C CP 2 : y 2 = z 3 + f(x)z + g(x),x G CP 1 , 

and 

S v C CP 2 : y 2 = x 3 + f(z)x + g(z), z G CP 1 . 

Here z, y are nonhomogeneous coordinates in CP 2 , with canonical symplectic 
form f2 = dzhdx and similarly x, y are nonhomogeneous coordinates in CP 2 , 
with the canonical symplectic form Q = dx/ ^ dz . 

Now the local action coordinate I{x) is computed |J as 

dl(x) = T(x)dx = ^1/2ti <j> dz/y^j dx, 

where [A x ] G Hi{£ X) 1) is a chosen A-cycle. Analogously, the dual action 
coordinate I D (z) satisfies 



dI D (z) = T D {z)dz 



1/27T j) dx/y^j dz, 
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where [L z ] G #i(C 2 ,Z). 

Now we have the following chain of transformations: 

dz Adx dx A dz 

y y 

dz dl dx dI D 
A -— = A 



V T(x) y T D (z)' 
or 

dz dx 

A dl = —p- ; A dJ 



2/T(x) 

Using the expressions for the angle variables d(p = y ^ z ^ and d(p D = yT d D^ 
we obtain 

dcp A dl = -dip D A dI D . 
Keeping in mind the interpretation of ||, the equality above amounts to 
K{k)dp Jac A dl = -K{k)dp JacV A dI D , 

or 

dP AdQ = -dp A dq, 

the starting point of f|. 

Remark The duality transformation x — > z, z — > a; between two elliptic 
pencils can naively be interpreted as a local "mirror transformation" between 
two "1/2K3" surfaces PS| : thinking of (z,x) of the first cubic as the Kahler 
and complex parameters we, obtain that (x, z) on the second becomes the 
corresponding pair of parameters. Hence the duality at hand changes the 
"Kahler" and "complex" structures of the cubic surfaces. 

3.4 Two coupling constants 

There exist natural generalizations of the double-elliptic system so far de- 
scribed involving additional parameters expressing noncommutativity. In 
terms of the intersection of quadrics the number of the independent coupling 
constants amounts from the number of their independent coefficients. One 
simple case involves such counting for the noncommutative T 4 . In this case 
all couplings follow from noncommutativity, hence the maximal number of 
parameters can be counted from the antisymmetric matrix dij defining the 
Poisson bracket 

{xii Xj} Oij (27) 

and equals to six In the stringy framework the couplings correspond 

to the fluxes of rank-two field switched on along different directions in four- 
dimensional manifold. 

One interesting example of two coupling constants comes from the fol- 
lowing system of quadrics 

S~l 2 ,y. 2 I ,-yj2 , y» 2 "1 

y Ju ^ \^ ju ^ J. 2 — J- 

g 2 x\ + x\ — x\ = k 2 
—g 2 x\ + x\ + x\ = 1 
— g 2 x\ + x\ + x\ = k~ 2 
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which is written in the most symmetric form. This system yields the following 
set of the Poisson brackets 

{xi,x 2 } = g 2 x 4 x 5 x 3 x 6 , {xi,x 3 } = g 2 x 4 x 5 x 2 x 6 , {x 2 , x 3 } = 

{x 5 ,Xi} = -X 2 X 3 X 4 X 6 , {X 5 ,X 2 } = -2i2 3 2 4 2 6 , {X 5 ,X 3 } = -£1X224X6 

{25,24} = -g 2 x 1 x 2 x 3 x 6 , {2 5 ,2 6 } = 0, {x 6 ,X4} = -g 2 x 1 x 2 x 3 x 5 . 

Now if 25 is taken as the Hamiltonian then the coupling g corresponds to 
the "interaction" of the coordinates, while g corresponds to the "interaction" 
of momenta. The duality now acquires a more symmetric form and corre- 
spond just to the simultaneous interchange of x 2 and 25 and the couplings. 
One could also consider examples of self-dual systems by choosing, say, 2225 
or x 2 + 25 as Hamiltonians. 

4 Degenerations 

Let us briefly consider possible degenerations and limits of the model. If one 
sends k to zero the dynamical system degenerates to the elliptic Ruijsenaars 
model and the corresponding four dimensional manifold reduces to [T 2 x C*] 9 
which will be explicitly described below. At the next step one can reduce the 
model to the elliptic Calogero model with the manifold [T 2 xC] 9 . Alternately, 
if k is sent to zero the system reduces to the one which is dual to the elliptic 
Ruijsenaars and Calogero models respectively. 

The constructions involving quadrics provide a very explicit description 
of the noncommutative manifolds above. Indeed let us consider limit k — > 
which amounts to the system of three quadrics in C 5 

1 2 — 

,71 2 , yi 2 — ^,2 

1 3 — 

-# 2 2 2 + 2 2 + 2 2 = 1. 

Geometrically this looks like a cylinder which fibered over the elliptic curve. 
The corresponding deformed four- manifold for the Calogero system [T 2 xC] 9 
is a "pinched" CP 1 fibering over the elliptic curve. 

Let us discuss the relation with the Sklyanin algebra 0. To this end let 
us interpret the Poisson bracket relations in the Sklyanin algebra 

{S a , S } = 2{Jp — J 7 ) SpS^ (28) 

{S a , Sp} — 2SoSj (29) 

as an example of the Mukai-Sklyanin algebra generated by two quadrics in 
CP 4 . The quadrics 

2Qi = ELi s 2 n (30) 
2Q 2 = S 2 + ELi JnS 2 n (31) 

coincide with the center of the Poisson bracket algebra. Hence the Sklyanin 
algebra fits into the general scheme. 
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Now recall the observation made in |30| that the Hamiltonian of the ellip- 
tic Ruijsenaars system coincides with the generator of the Sklyanin algebra 
So. On the other hand, following our approach, the double-elliptic Hamil- 
tonian coincides with the coordinate x$. This means that in the limit of 
vanishing k the system of four quadrics in CP 6 reduces to a system of two 
quadrics in CP 3 corresponding to the Sklyanin algebra. 

The degeneration to the trigonometric Ruijsenaars model can be per- 
formed in a similar way and the corresponding Hamiltonian can be expressed 
in terms of the coordinates which follow from the realisation of U q (SL(2)) in 
terms of the intersection of two quadrics in C 4 . The explicit degeneration in 
this direction goes as follows. Let us first consider the algebraic structure in- 
termediate between the Sklyanin algebra and U q (SL(2)). Following |TL|, we 
define generators A, B, C, D of the degenerate Sklyanin algebra at "small" 
h = e i7rT (^r -> 0): 

h~ 1/2 

A = - — — - — (cos tit] S -Msin 7rr]S 3 ), 
2 sm Z7r?7 

D = - — — - — (cos 7ii] S - i sin irr] S 3 ), 
2 sm 2-kt] 

C = (S 1 -iS 2 ), 

2 sin 2ttt] K lh 

h -3/2 

B = -——(S 1 + iS 2 ). (32) 
8 sm 2Tirj 



The generators A, P, C, D satisfy the quadratic algebra [3jJ : 

DC = e 2mv CD, CA = e 2niv AC, 
AD- DA = -2i sin 3 2-kt] C 2 , 
A 2 -D 2 

BC - CB 



2% sin 2717? 

AB - e 2wir 'BA = e 2mn DB — BD = % - sin 4^77 [CA - DC). (33) 

The Casimir elements are the following quadrics 
Q 1 = e ^AD -sin 2 2tt77 C 2 , 



Qi = L.. PC -cos 2nr]C 2 . (34) 



e -27rir/^2 _|_ ^iriri Jj2 

4 sin 2 2txt\ 

and this algebra fits within the general formula for a generalised Mukai- 
Sklyanin algebra. Some generalisations of this picture for the trigonometric 
many-body system can be found in [|32|. It would be extremely interesting 
to present the elliptic Ruijsenaars-Schneider many-body systems explicitly 
in a similar manner in terms of the quadratic elliptic algebras of |33| or their 
tensor products. 

To get the U q (SL{2)) itself from the algebra above it is necessary to 
consider the contraction e — > 0: 

P -> P, A -> eA, D — > eD, C -> e 2 C. 
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and take the following quadrics given by the Casimirs 



Qi = AD 

-2/12 , 2 f)2 

where q = expinrj. The Poisson brackets following from the algebra of the 
quadrics reproduce the standard Poisson structure of U q (SL(2)). 

One further remark concerns the degeneration to the Calogero model. 
In principle we could obtain this by two paths. First, we can represent the 
system in terms of the quadratic Hamiltonian with linear Poisson bracket. To 
do this it is useful to consider the limit of large value of the Casimir Q 2 and 
express So in terms of other generators of the Sklyanin algebra. Expanding 
the root to first order we immediately get the representation of the elliptic 
Calogero system as an elliptic rotator 

3 

H = J2 JiSf (36) 

n=l 

which reduces to the standard form after a redefinition of the coordinate. 
However we could also consider the quadratic algebra using the degeneration 
of the action of x on the function f(u) 

f( v 6>iifa)(M2tt ~ ly)f(u + y)- 9 n (-2u - l)f(u - r,)) 
Xof(u) = W (37) 

where / is the spin of representation corresponding to the coupling constant 
in the dynamical system. To get the elliptic Calogero system we have to send 
77 — > hence the difference operator becomes the differential one. 

For completeness let us note that the periodic Toda two-body system can 
also be described in terms of a top on the quadratic algebra. To this end 
consider the following Poisson algebra 

{A, A 2 } = 2A 3 (4 - A 2 ), {A 3 , A 2 } = A 2 , {A 1} A 3 } = A l + Aj. 

Then the Hamiltonian of the Toda system is the linear function 

H Toda = l/2(A 1 + A 2 ) (38) 

This picture can also be considered as the further degeneration from the 
elliptic Calogero model via the Inozemtzev limit. 



5 Other incarnations of the double-elliptic sys- 
tem 

5.1 Double-elliptic systems and decoupled Nahm tops 

We are now able to give another polynomial description of the double-elliptic 
system observing that it has the form of Fairlie's "elegant" integrable system 
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35fl for n = 4. Ignoring the unimportant coordinate Xq we have, 



x\ = X2X3X4 (39) 

x 2 = £1X3X4 (40) 

x 3 = X1X2X4 (41) 

±4, = g 2 x 1 x 2 x 3 . (42) 

This system admits a beautiful description as a decoupled system of Euler- 
Nahm tops after the following change of variables: 

u + = X3X4 + gxix 2 , v + = x 2 x 4 + gx x x 3 , w + = x x x 4 + gx 3 x 2 
u_ = X3X4 — gx\x 2 , v_ = X2X4 — gx\x 3) w_ = x\x 4 — gx 3 x 2 - 

In the terms of the new variables the double-elliptic system (|39| ) is equivalent 
to 

u + = v + w + (43) 

i) + = w + u + (44) 

w + = u + v + . (45) 

and to 

m_ = v_W- (46) 

?>_ = W-U~ (47) 

w_ = U-V-. (48) 

Geometrically this change of variables means a passage from the intersection 
of the quadrics to the direct ("decoupled") product of two elliptic curves 
£ + x £_ given by the Casimirs of the models ([l"3 ) and 



S + :u 2 + -v 2 + = k 2 (E 2 -l) 
u\-w\ = (k 2 -2)(E 2 -l) 

and 

E_ : u 2 _ - v 2 _ = k 2 (E 2 - 1) 
u 2 _-w 2 _ = (k 2 -2)(E 2 -l). 

This result is reminiscent of the result of of Ward |M] that the second 
order differential operator with Lame potential n(n + l)/2sn 2 (q\k) can be 
factorized into a product of the first order matrix operators (d + A)(d — A) 
with the matrix A = (Ai, A 2 ,A 3 ) satisfying the Nahm equations 

A-i = Gijk [Aj , Ak\ . 

This property manifests a hidden supersymmetry underlying the quantum 
mechanics with Lame potential. The factorization of the double-elliptic sys- 
tem into two decoupled tops suggests the existence of a hidden SUSY in this 
case too. 
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5.2 The Double-elliptic system as an example of a Nambu- 
Hamilton system 

Our polynomial parametrization of the double-elliptic system also provides 
and example of a Nambu-Poisson structure. Such a structure was introduced 
by Nambu in 1973 |38] as a natural generalisation of the Poisson structure 
(which corresponds below to n — 2). Such structures were recently exten- 
sively studied in ||39| . Recall that a Nambu-Poisson manifold M n is is a 
manifold endowed with an antisymmetric n- vector rj 6 A n (TM). The field 7] 
defines an n-ary operation: 

{,...,}: C°°(M) m i — ► C°°(M) 

such that the three properties are valid: 

1. antisymmetry: 

{/l, fn} = (-1) CT {/<t(1), fa(n)}, ® € S n ] 

2. coordinate- wise "Leibnitz rule" for any h 6 C°°(M): 

{flh, f n } = fl{K f n } + h{f t , fn}] 

3. The "Fundamental Identity" (which replaces the Jacobi Identity): 

{{fl, fn}, fn+l, fln-l} + {fn, {fl, {fnY fn+l}, fn+2, /2n-l}+ 
+ {fn, •••) f2n-2, {fl, fn-1, /W-l}} = {/l; fn-1, {fn, /Wi-l}} 

for any A,.., / 2n _! eC°°(M). 

Dynamics on a Nambu-Poisson manifold is defined by n — 1 Hamiltonians 
...,H n _i giving the Nambu-Hamiltonian system 

-77 = {H\,H 2 , ...,H n -i,x}. 

The most common example of a Nambu-Poisson structure is the so-called 
"canonical" Nambu-Poisson structure on C n with coordinates Xi, ...,x n : 

{fi, ■-, fn} = Jac{f x , f n ) = "•'■/") 

d{x 1 , ...,x n ) 

Our constructions involving quadrics give non-trivial examples of a Nambu- 
Hamilton dynamical system. Namely, the system of three quadrics in C 5 

T2 , yi 2 ~1 ,-y* 2 ,yi 2 j^*2 j-j 2 ™ 2 I ™ 2 I ,->i 2 "1 

ju <2 — J- ^ *Ju -y 3 — ? U 1 ~~ 4 r" 5 — 

admits a section by the choice of the level x^ = E and the 1-parameter 
intersection of three quadrics in C 4 

Q 2 = - Xg = /c 2 

define a Nambu-Hamilton system, which is our double-elliptic system: 

—jt = {Qi^Q2,Qz,Xi}. 
at 
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6 N-body systems 



In this section we discuss the challenging problem of constructing the many- 
body generalisations of double-elliptic systems. We will not present an ex- 
plicit realisation of such systems, but some problems and proposals concern- 
ing the matter will be given. Let us emphasise that any finite dimensional 
two-body system can be generalised in at least in two directions. 

The first one involves spin chain type systems. The starting point to be 
generalised is the Sklyanin Lax operator expressed in terms of the generators 
of the Sklyanin algebra obeying the RLL = LLR relations with Belavin's 
elliptic i?-matrix. The corresponding iV-body system is the XYZ spin chain. 
Now we would like to construct the local Lax operator expressed in terms 
of the generators of a generalised Mukai- Sklyanin algebra of n quadrics for 
n > 2. It can be shown that it is impossible to introduce a local Lax operator 
whose matrix elements are such generators assuming that Belavin's R matrix 
intertwines it. Therefore to pursue this line of generalisation one has to 
develop some generalisation of the elliptic i?-matrix. 

Another line of generalisation would be to seek the putative many-body 
double-elliptic system as a generalised Hitchin system. Hitherto the most 
complicated system treated in this way has been the elliptic Ruijsenaars 
model which can be obtained via the Hamiltonian reduction procedure from 
the following moment map equation |16| 



g(z)h(z)g-\z)h- 1 (qz)=fji (49) 

where g and h are elements of affine SL(N) and /i is the fixed level of the 
moment map. 

To obtain the double-elliptic system one presumably has to solve the 
following equation 

g(z)h(z)g- 1 (qz)h- 1 (qz)=ti (50) 

which can be considered as the moment map equation for the Heisenberg dou- 
ble of the double affine algebra. The symplectic manifold to be constructed 
is the generalisation of the cotangent bundle to the double affine group which 
provides the phase space for the elliptic Ruijsenaars model. If the solution 
to this equation could be found (in some particular gauge) the Hamiltonian 
of the many-body double-elliptic system would be given by J Trg(z). The 
results obtained in |32| are of help for this approach. 



To link the Heisenberg double approach with the geometrical consider- 
ations of this paper it is instructive to adopt the fundamental group inter- 
pretation of the moment map equations. The simplest equation of such type 
corresponding to the trigonometric Ruijsenaars model looks like 

ghg^h- 1 = fi. (51) 

This has no coordinate dependence and just represents the relation following 
from the fundamental group on the torus with one marked point. To discuss 
the double-elliptic case in a similar manner we have to consider connections 
on the whole four dimensional manifold. Actually we are interested in a 
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twisted bundle which in the simplest example of the torus results in t'Hooft's 
consistency condition for the twists 



VLi(z + a,i)Vtj(z + aj)Qi l (z + cn)Qj 1 (z + a,) = Z t j. (52) 

Here Zij belongs to the center of the group and represents the matrix of 
fluxes. The torus case is too simple to generate the double-elliptic system. 
However if we solve (|5^) then the hamiltonian j TrQi would provide a very 
degenerate example of such a "double-elliptic " system. 

One more powerful approach is based on the separation of variables pro- 
cedure. Since we have identified the phase space of the two-body double- 
elliptic system one could imagine that the Hilbert scheme of the n-points on 
this manifold would solve the problem. However, a serious technical difficulty 
arises since there is no effective description of Hilbert schemes of points on 
generic four dimensional manifolds. A related issue involves the identification 
of integrable many-body systems in terms of noncommutative instantons El 



To obtain the description of the double-elliptic system in terms of noncom- 
mutative instantons the ADHM description of instantons on commutative 



and noncommutative K3 manifolds should be developed. (See |42| for steps 
in this direction.) 



7 The relation to D=6 SUSY gauge theories 
with massive adjoint matter 

Let us briefly discuss the possible application of double-elliptic systems. It is 
known that there is the map between classical integrable many-body systems 
of Hitchin type and their generalisations (|T5|, 18 1 and the references therein) 
with supersymmetric gauge theories in d = 4, 5, 6 dimensions with the dif- 
ferent matter content. This relation can be briefly formulated as follows. 
The solution to the classical equation of motion of the integrable system 
provides the spectral curve embedded into the higher dimensional compacti- 
fication manifold. The gauge theory is defined on the worldvolume of the M5 
brane wrapped around the spectral curve. The particular integrable system 
is in one to one correspondence with known supersymmetric gauge theories. 
Namely, the space of integrals of motion is mapped into the Coulomb branch 
of the moduli space in the gauge theory, while the modular parameters of 
the spectral curve of the integrable system amount to the running coupling 
constants on the field theory side. 

Previously the elliptic Ruijsenaars integrable system, of Hitchin type, was 
identified with an N = 1 d = 5 SUSY gauge theory with one compact dimen- 
sion and adjoint hypermultiplet pi 0, [2"1]|. The double-elliptic systems we 



are looking for presumably correspond to the next step, namely an iV = 1, 
d = 6 gauge theory with adjoint hypermultiplet when two dimensions are 
compactified. There are two natural elliptic moduli in the gauge theory: one 
is the complexified coupling constant while the second is the modulus of the 
compact torus in the fifth and sixth dimensions. These are to be identified 



23 



with the moduli k, and k in the quadrics above. The coupling constant g 
as before can be related to the mass of the adjoint field. We suggest that 
the genus 2 curve (^) we have found for the two-body double-elliptic sys- 
tem plays the role of the Seiberg-Witten curve providing the renormalised 
coupling constant in d = 6 gauge theory with two compact dimensions. 

The Coulomb branch of the moduli space of this theory can be also de- 
scribed, via an additional compactification of one further dimension, and 
consideration of the resulting 2 + 1 theory with N = 4 SUSY as the theory of 
NS5 branes compactified on T 3 |^] . The moduli space of this theory can be 
mapped into the moduli space of noncommutative instantons on T 4 and K3 . 
The noncommutativity corresponds to the fluxes of B field along some direc- 
tions. Generically fluxes can be switched on in all directions, corresponding 



to the complete set of twists in the language of [}12j1 . The system considered 
in this paper corresponds to a flux only along the fiber of the elliptic fibration 
which plays the role of the coupling, constant in the dynamical system and 
the mass of adjoint in the field theory. To link this with the integrable system 
one could adopt the description of the instanton moduli space via the spectral 
curve endowed with the spectral bundle [i^j. This description emerges after 
a T duality transformation of the initial four dimensional manifold along two 
dimensions. 

It has been argued that the instanton spectral curve coincides exactly 
with the Seiberg-Witten curve describing the Coulomb branch of the moduli 
space of some N = 2 theory. The reason for this relation is that moduli 
space of the instantons on noncommutative T 4 is equivalent to the moduli 
space of the little string theory on T 3 . The twists correspond to the structure 
of the T 2 fibration over the the base T 2 . In the generic situation there are 
twists over the fiber and base tori. Geometrically the situation matches with 
the picture of the intersection of the quadrics above. The coupling constants 
correspond to the twists and it is therefore natural to conjecture that the 
double-elliptic n body system (or possibly a degenerated version) is related 
to n, U(l) instantons on noncommutative T 4 . 



8 Conclusion 

In this paper we have elaborated in some detail the simplest double-elliptic 
systems and their relation to the generalised Mukai-Sklyanin finite dimen- 
sional algebras. The previously discovered Hamiltonians for these systems 
were placed in an algebro-geometric setting allowing various connections with 
Nahm tops to be made. 

The generalisation to the many-body case still remains a challenging prob- 
lem. We have however identified several mathematical problems to be solved 
to make progress in this direction. Namely, it is highly desirable to find the 
generalisation of the i?-matrix related to the generic algebra of quadrics, to 
develop the ADHM description of Hilbert scheme of points on intersection 
of the quadrics and to study further the ADHM description of instantons on 
noncommutative K3 manifolds. It is also necessary to find the U SU(N)" gen- 
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eralisation of the algebra of quadrics which would involve a larger number of 
generators but keeping the power of the generators involved (cubic, quartic 
e.t.c) the same as for "SU(2)" case. Such a generalisation is known |33| only 
for the Sklyanin algebras of 2 quadrics. The corresponding quadratic Feigin- 
Odesskii algebras involve an arbitrary number of generators and can be built 
from the vector bundles on the elliptic curve. The last problem to be men- 
tioned is that of finding a precise formulation of the notion of the Heisenberg 
double to a double affine algebra. It would be also very interesting to com- 
bine the ideas from this paper with the recent two-dimensional generalisation 
of the Calogero and Hitchin type systems suggested in H44|,[fl5[ . 

Besides the problems above, which can be considered as primarily math- 
ematical ones, there are a number of issues relevant for the physics of D=6 
gauge theories. First, let us emphasise that the genus two spectral curve fllBD 
provides the low energy effective action and the spectrum of BPS particles in 
D=6 theory. It would be interesting to pursue further the duality arguments 
relating the theories in different dimensions. One more interesting possibility 
along this way involves the idea of the deconstruction of dimensions, starting 
from D=4 theory. The deconstruction of D=6 theories with two compact 
dimensions has been considered recently We expect that the double 

elliptic system considered here could be derived from the "double lattice" of 
the integrable models corresponding to D=4 if the nonperturbative effects 
would be taken into account. This correspondence has been proven for the 
D=5 models with one compact dimension which can be derived from the 
"lattice" provided by the quiver type D=4 theory p| . 
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10 Appendix 

We give for completeness the computations showing that the double-elliptic 
Hamiltonian system in canonical variables is rewritten in the form (|2|) under 
the change of variables (3-8). 
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10.1 A Hamiltonian equations 

We begin with the following relations 
dsn(u\k) cn(u\k) dn(u\k) 



dk kk' 2 

dcn(u\k) sn(u\k) dn(u\k) 

dk kk' 2 
ddn(u\k) k sa(u\k) cn(u\k) 



[-E(u) + k' 2 u + k 2 su{u\k) cd(u\k)] 



[E(u) - k' 2 u - k 2 sn{u\k) cd{u\k)] 
[E(u) - k' 2 u - sn(u\k) dc{u\k)] 



dk k' 2 

where k' 2 = 1 — k 2 , with cd(u\k) = and so forth. Here E{u) = E(u, k) 

is the elliptic integral of the second kind 

/sn(it|fc) /]_ _ ^2^2 fu 
yj — — p-dt = J dn 2 (v\k)dv. 

Only the first of these relations needs to be established, with the latter 
two following from partial differentiation of the identities 

sn 2 (w|£;) + cn 2 (u\k) = 1 

and 

dn 2 (u\k) + k 2 sn 2 (u\k) = 1 

respectively. 

To establish the first we take the partial derivative with respect to k of 
the definition 

r sn(«|fc) dt 



u 



^/(l-t 2 )(l-k 2 t 2 ) 



giving 



1 dsn{u\k) f sn{ulk) t 2 dt 



cn(u\k) dn(u\k) dk J (1 - kH 2 )y/(l - t 2 )(l - k 2 t 2 ) 

Upon change of variable we have 

1 dsn{u\k) f u sn 2 {u\k) 

~ (JjUj 



k cn(u\k) dn(u\k) dk J dn 2 (u\k) 

which is straightforwardly integrated to yield the first identity. 
Let us set 



g 2 _ / g 2 k 2 a _ 

" = Xl sn 2 (q\kY P ~ V sn^(#)' /T 

The equations of motion for the Hamiltonian H — a(q) cn(/3(q)p\A) yields 

dH 

q = — — = -a(3 sn((3p\A) dn(/3p|A) 
dp 
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and 

OH 

p = — — = —d q a ca(/3p\A) + pa sn(f3p\A) dn((3p\A)d q (3 — ad A cn(f3p\A)d q A. 
Here 

d A - -a 2 ~k(l - p ^n(g|A;) dn(g|fc) 

We note that p is invariant under p — > p + even though the explicit 
momentum terms and the elliptic integral of the second kind in the last term 
individually are not. 

Now introduce the quadrics 

1 2 — 

/ y» — ,->-• 2 7^2 

1 3 — 

a 2 = 1 + g 2 a: 2 

which may be realised by 

1 cia(q\k) dn(q\k) 

1 sn(g|A;)' 2 sn(g|A;)' 3 sn(g|A;) 

Consider further the quadrics 

vl + yl =i 

realised by 

?/i = sn(/3p|A), y 2 = cn(/3p|A), y 3 = dn((3p\A). 
Upon setting 

x 4 = ayi, x 5 = ay 2 , x 6 = f3y 3 
we recover (^j and find that H — x%. Further 

x 4 = sn(f3p\A)d q aq + a cn(/3p| A) dn(/3p| A) (pd q f3q + + «(9a sn(/?p| A)<9gAg. 
We see 

x 4 = (yf 2 xix 2 a;3a; 6 

For this we don't actually need the explicit expressions for d A sn((3p\A) and 
c?A.cn(/3p|A) for these appear in the combination 

-a 2 /3<9 9 Adn(/3p|A) (sn(/3p|A)<9 A sn(/3p|A) + cn(/3p\A)d A cn(f3p\A)) 

which vanishes as a result of the first identity given earlier. 
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10.2 B Canonical 2-form in different parametrizations 

Let us compute the expressions for the canonical 2-form dp A dq given in the 
local coordinates of the initial elliptic curves as a 2-form on the phase surface 
in the different parametrizations. 

We begin with the following identity in the region x 5 7^ 

_ , , dx\ A dx^ 
\l = dp A dq = . 

x 2 x 3 x 5 x 6 

Indeed, we have from the obvious relations 

dH A dq dx^ A dx\ 
dp Adq = — 7T77 — = . 

^ x 2 x 3 x i x 6 

and (f|),( |IO|) and (|i~l~D amount to dx 5 A dx\ = ^dxi A dx A . 
In the notation of Section 2 

dX Adz 

This follows from the identities 

dzf\dx\ dzAdx dzAdx dzAd(l/X) dXAdz 
x 2 x 3 x A x 6 x l x 2 x 3 x 4 x 6 yx 6 Y/X 2 x 6 Yx 6 
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